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Transparency of the metal/superconductor interface is discussed in a frame of the Kupryanov- 
Lukichev (K-L) quasiclassical boundary conditions. It is shown that the original K-L boundary 
transparency is expressed via the single-channel Landauer conductance. It is a source of dramatic 
discrepancy between the theory and experiments on the superconducting proximity effect. The 
improved multi-channel derivation is proposed for the boundary conditions, which rehabilitates con- 
sistency of the quasiclassical approach to the boundary transparency problem, and eliminates at 
least an order of magnitude of discrepancy between the theory and the experiment. The influence 
of the interface roughness and wave function symmetry mismatch at the boundary on the inter- 
face transparency is discussed. Warning is made to theories that use the quasiclassical boundary 
conditions with the interface transparencies close to unity. 



Since works by Eilcnberger Larkin and Ovchin- 
nikov and Usadel Q the quasiclassical theory is 
main calculation approach in physics of superconduc- 
tivity. Most of calculations utilize the Usadel theory 
with the Kupriyanov-Lukichev (K-L) boundary condi- 
tions because the short mean-free path approxima- 
tion matches most experiments with conventional super- 
conductors (note that the short coherence-length, high- 
temperature superconductors refer mainly to a clean- 
limit case). As long as the boundary resistance, which 
enters the K-L boundary conditions, is considered as a 
phenomenological fitting parameter, the Usadel theory 
gives adequate description of contact phenomena that 
has been demonstrated, for example, for the Josephson 
effect by Golubov et al [5j and for the proximity effect by 
Buzdin Ifl. However, it was pointed out by the Salerno 
group [21 |S| that problems immediately arise when one 
tries to reproduce the best-fit boundary resistance (or 
transparency of the interface) , obtained from proximity- 
effect experiments, making use of the physical parameters 
of materials in the contact. The authors employed the 
band offset (Fermi momenta mismatch) model for the 
quantum-mechanical transmission coefficient, 
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which is plausible approximation as far as the free- 
electron model is used as a background for the quasi- 
classical theory (vi x is a projection of the electron Fermi 
velocity of the i-th metal on the direction perpendicular 
to the interface) . The quantum-mechanical transmission 
coefficient D is related to the boundary resistance Rb via 
parameter 7b as follows 
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where A is an area of the contact, a 2 = e 2 p 2 l 2 /3ir 2 ti i 

1 /2 

is the normal-phase conductivity, £2 = (1iD 2 /2ttT c ) is 



the coherence length, l 2 and p 2 are the mean free path 
and the Fermi momentum of electrons, D 2 = v 2 l 2 /i is 
the conduction electron diffusion coefficient, T c is the su- 
perconducting transition temperature, and, finally, x 2 — 
cos #2, where 9 2 is the quasiparticle incidence angle. The 
subscript "2" refers to a metal #2. The parameter 75 
varies between zero (virtual interface between identical 
metals) and physical infinity (an ultimate case of highly 
resistive insulating barrier at the interface). Following 
Aarts et al. @ the Salerno group characterized the in- 
terface by the transparency parameter, 
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which varies between and 1. From the fitting of their 
data on the proximity effect in Nb/Pd system the 
authors obtained 7^a;p(Nb/Pd)~0.46, whereas their esti- 
mation from Eqs. Q and © gives 7^ eor (Nb/Pd)~0.98 
using the Fermi velocities 1^=2.73 xlO 7 cm/s and 
w Pd =2.00xl0 7 cm/s. For Nb/Cu and Nb/Ag cou- 
ples the Salerno group obtained T exp (Nb/Cu)c^0.30 
as compared with 7^ eor (Nb/Cu)~0.5 (wcti = l-57xl0 8 
cm/s); and 7^, xp (Nb/Ag)~0.33 as compared with 
% eor (Nb/Ag)~0.55 (im 9 =1.39x10 8 cm/s) 0. Similar 
low interface transparencies have been derived from prox- 
imity experiments on Nb/Al couples: 7^ (Nb/Al) ~0. 2- 
0.25 (Ttheor -0.8) and Nb/CuNi couple: 



T e:C p(Nb/CuNi)~0.2-0.35 



One may derive from 



the above data that there is systematic and drastic dis- 
crepancy between the experimentally established trans- 
parency of the S/N interface and the quasiclassical the- 
ory which is applied to reproduce the transparency data. 
The discrepancy looks even more dramatic if we recognize 
that one needs unrealistic difference in the free-electron 
Fermi energies, one hundred times and more, to repro- 
duce the boundary resistance, say, Nb/Cu couple. The 
interface roughness or mismatch of the wave functions 
symmetry at the interface can not compensate so dras- 
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tic misfit between the experiment and the theory. Then, 
one may come to a conclusion that there is global incon- 
sistency in the quasiclassical approach when applied to 
description of the contact phenomena. 

In this paper we show that source of the discrepancy 
is a very approximate treatment of the boundary resis- 
tance, which has been made upon derivation of the K-L 
boundary conditions. We develop further the K-L ap- 
proach and propose new derivation of the boundary con- 
ditions for the Usadel equations which relaxes mainly the 
above mentioned discrepancy. We argue that approxima- 
tion of highly transparent interface (the quasiclassic su- 
perconducting pairing function is continuous across the 
interface), which was oftenly used in calculations of the 
contact phenomena, is hardly (if ever possible) to realize 
in an experiment with conventional superconductors. 

We start from the Zaitsev boundary conditions 
for the Eilenberger equations following Kupriyanov and 
Lukichev: 

9al = 9a2 = 9a, (4) 



(cos 9 f a ) =const(z), obtained after the solid-angle aver- 
aging of equation JSJ, and the antisymmetric function 
f a2 ~ — l 2 cos 9 2 (dF 2 / dz) , found from Eq. J7J) at z 3> 1% 
in the lowest order on anisotropy, we obtain from Eq. Q 
the first linearized boundary condition at z = 0, 

, 2^Fi 2 dF 2 . . 
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where Fi — (/j) is the isotropic anomalous Usadel func- 
tion. 

To derive the second boundary condition we divide at 
first both sides of equation © by 2R: 

fa = ^ (/.a - fsi) ■ (10) 

Then, we substitute the approximate antisymmetric 
function, utilized above, for the left-hand side, replace 
the symmetric Eilenberger functions in the right-hand 
side by their Usadel counterparts, f s i — ► F si , multiply 
HI U| l by cos 9 2 and average over the angle theta. This pro- 
cedure yields the second linearized boundary condition, 



g a ji? (1 - g a g a ) + — (g sl - g s2 ) j 

= j (9si ~ 9s2) (9si + 9s2) , (5) 

where g s ^ a )i are the symmetric (antisymmetric) Green 
functions, and R = 1— D is the reflection coefficient. The 
boundary condition J3J ensures continuity of the charge 
current at the interface, while the second one, Eq. JSJ, 
relates the quasiparticles flux to a drop of their density 
at the interface, and depends on transparency of the in- 
terface. Using the explicit matrix structure of the Eilen- 
berger functions |l3j close to the transition temperature 
T c , we arrive at the basic for our discussion, linearized 
Zaitsev boundary condition 

2Rf a = D(f s2 -f sl ), (6) 

where fi is the anomalous pairing Green function, which 
describes superconductivity in the system. We have to no- 
tice here that equation © has been derived without any 
approximation concerning transparency of the interface. 
To proceed further we need also the linearized Eilen- 
berger equations for the normal metal which read: 

i cos e 2 ^i + f a2 = o, (7) 



lcos9 2 ^+f s2 = (f s2 ), (8) 
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where the angular brackets mean averaging over the solid 
angle. The original logics of K-L derivation is as fol- 
lows (see also review 0], Ch. 4): using a constraint, 



676-^=^1-^2, (11) 

with 7b given by Eq. J3J . Equations 10 and are well 
known linearized Kupriyanov-Lukichev boundary condi- 
tions Q for the Usadel functions. It has been shown by 
Lambert et al 0] that they are valid only at low trans- 
parencies of the S/N interface. Why does it so despite of 
the generic Zaitsev BC © are valid for arbitrary trans- 
parency? 

To answer the question let us have a look at the phys- 
ical meaning of the linearized Zaitsev BC (|10f> which is 
parental for the K-L one, Eq. It has a sense of 

the boundary condition for a particular trajectory deter- 
mined by the angles <p and 9. The language of trajectories 
is suitable when there is a continuous domain of quasipar- 
ticle incidence angles which accommodates macroscopic 
number of trajectories. In a case of a very few travers- 
ing trajectories it is more suitable to use a language of 
quantum- mechanical conductance channels. Then, the 
single- channel Landauer conductance |lfil | (in units of the 
conductance quantum e 2 /irh) can be immediately recog- 
nized in the ratio D /R which enters the right-hand side of 
equation (|10() . Subsequent derivation procedure involves 
integration over the incidence angles of the single-channel 
Landauer conductance. Implication of the single-channel 
conductance to BC is crucial for understanding the source 
of insufficiency of the original K-L approximation repro- 
duced above. 

Indeed, the integrated single-channel Landauer con- 
ductance matches physical expectation, that conduc- 
tance of the interface is zero when the transparency coef- 
ficient D is zero (boundary resistivity is infinite). Alter- 
natively, the conductance is infinite (boundary resistivity 
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is zero) when the contacting materials are identical, and 
the interface is virtual. On the other hand, practically 
in all experiments with superconducting heterocontacts 
we have many channels for conductance through the in- 
terface. From the physical point of view one may expect 
that the multi-channel conductance should enter the cor- 
rect formulation of the boundary condition Then, 
we have to answer on a next question: is the angle- 
integrated, single-channel Landauer conductance repre- 
sent correctly the multi-channel conductance? The an- 
swer "no" is given at the end of the Section II of the paper 
by Buttiker, Imry, Landauer and Pinhas [TtI ] (hereafter 
will be referred as BILP). In the Section IV of their paper 
the authors derive the multi-channel conductivity of the 
interface, valid at arbitrary transparency: 
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and Vi are velocities of quasiparticles moving to the left 
(I) and to the right (r) over N and N' channels of con- 
ductance, respectively. For a flat interface of macro- 
scopic area the summations over the conductance chan- 
nels can be replaced by integration over the incidence 

angles (lif: Ei=\ ~~ * / dtp a dO a sin 9 a , and the integra- 
ls 

tion runs over the semi-sphere (SS) in the direction of the 
particles velocity. Dependence of the multi-channel con- 
ductance i|12|) on the Fermi-energies of contacting metals 
and roughness of the interface between them has been 
analyzed in detail by Garcia and Stoll [l8j |. 

To make evident relevance of a sheet conductance arj 
(conductance per unit square) of an interface for formu- 
lation of BC we rewrite Eqs. © and (|llfl in terms of 
conserving fluxes: 
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is nothing else but the sheet conductance of the inter- 
face in the Kupriyanov-Lukichcv approximation. How- 
ever, as we deduced from the above analysis, the ulti- 
mate formulation of BC (|15fl has to be expressed via the 



multi-channel sheet conductance (at the current status, 
with the last parentheses in Eq. I|16|) replaced by the 
parentheses from Eq. I112[0 . Now, the main direction of 
logics is to find an analogue or a good approximation of 
the multichannel conductance remaining within the qua- 
siclassical formalism. 

As the first step let us look what will happen if we 
withdraw transformation of the original, particular tra- 
jectory Zaitsev BC © into the single-channel BC form 
1)10(1. Starting from Eq. JBJ we replace: (1) the anti- 
symmetric Eilenberger function f a in the left-hand side 
by the approximate expression used upon derivation of 
Eq. @; (2) the symmetric Eilenberger functions in the 
right-hand side by their Usadel counterparts. Then, we 
have in an output 

dFo 

-2{1 - D)l 2 cos e 2 —^- = D{F 1 -F 2 ). (17) 

dz 

The procedure of the angular averaging described after 
Eq. HI U| l yields BC H15|) with the sheet conductance 
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The behavior of cr§ T1 on a ratio of the Fermi energies 
Rf = Sf2/sfi < 1 is given by the dash line in Fig. 1. It 
is easy to see by comparison with Eq. (|16f) that the main 
problem of K-L sheet conductance at high transparency 
of interface - the denominator, 1 — D, which is close to 
zero in the main domain of integration angles - is relaxed 
already in Eq. (|18f) . However, the first step does not take 
into account properly the influence of other conducting 
channels (or trajectories crossing the interface) on prop- 
agation of an electron along a particular trajectory. 

When there is a single conductance channel through 
the interface, all other bulk conducting channels ap- 
proaching the interface remain unperturbed. In our par- 
ticular case of the S/N junction the superconducting or- 
der parameter parameter is flat near the impenetrable 
interface with the only an atomic-size hole which acco- 
modates one conductance channel. In the case of the 
transparent interface, the other conductance channels in- 
duce mutual accomodation of the pairing function wave 
function and the order parameter from the both sides of 
the interface (see, for example, Fig. 3 of Ref. 0) creating 
gradients. Thus, the next step of approaching the true 
multichannel formulation could be consideration of the 
Eilenberger function gradients. 

Careful look at the derivation procedure shows that 
upon replacement of the symmetric Eilenberger functions 
in the right-hand side of © by their Usadel limits we 
have lost in Eq. (|17l) gradient terms of the same order 
that we have in the left-hand side of this equation. In 
the dirty limit, the electron mean free path is the short- 
est length compared with the superconducting or normal 
state coherence lengths. Then, in the spirit of Ref. [19j, 



4 



we expand the Eilenberger functions / s , in powers of the 
mean free path along the particular trajectory, keeping 
the gradient terms, 

f si — {f si} ~t~ Izi ~j — — Fi ~\~ ^zi , ; (^"^) 



and substitute them into the right-hand side of Eq. JHJ . 
After the angular averaging procedure we arrive at the 
second BC (tT3|) with the sheet conductance as follows: 




(20) 

Our sheet conductance (|2*0l is the next step to approx- 
imate the multichannel BILP sheet conductance l|12|) in 
the gradient approximation. The two terms in the de- 
nominator of the big parentheses in (|20|) are one-by-one 
counterparts of the similar terms, which we will have in 
(|12|l after replacement of Ri by 1 — Di in the denomina- 
tor. We believe (from considerations of energies involved 
in the transmission) that the result holds for the general, 
non-linear BC as well. Dependence of the sheet con- 
ductance l|2(JI) on Rp is displayed by the thick solid line 
in Fig. 1. It can be seen that the GT2 curve goes much 
closer to the BILP sheet conductance (dotted curve) than 
our first approximation given by Eq. I|18J) . The boundary 
conditions (|14|l and Ijl5(l can be expressed via the exper- 
imentally measured bulk resistivities of the contacting 
metals, p\ = a^ 1 and p2 = ovT 1 , an d the sheet resistance 
of the boundary, ARb = er^ 1 . 

We believe that our approach gives consistent quasi- 
classical approximation to the multi-channel BILP sheet 
conductance of the interface. The improvement is drastic 
in the region of high transparencies: our sheet conduc- 
tance ctq T2 is more than order of magnitude lower that 
the K-L one. Correspondingly, our sheet resistance of the 
interface, AR^ T = (cr§ T2 ) , in the region of R F < 1 
is more than order of magnitude higher than the K-L 
one. There are also other physical reasons that may in- 
fluence the boundary resistance. Garcia and Stoll have 
shown (Fig. 2 of Ref. 0) that the interface roughness 
increases boundary resistance, i.e. 7b . At the Fermi en- 
ergy ratio Rf = 0.4 this increase is about 20-60% for 
different models of the interface roughness. The reduced 
overlapping and symmetry mismatch of the wave func- 
tions may also essentially suppress transparency of the 
interface. Other physical aspect is that only the prop- 
agating (dispersive) part E prop of the total electron en- 
ergy, E to t — Ep r0 p + E rot , computed in the free electron 
model, contributes to the transport across the interface. 
So the narrow d-electron dispersion band may effectively 



produce reduced Rf ratios as soon as the localized ro- 
tational contribution E prop is subtracted from the total 
energy of the d-state. 

In a consequent way, one may think that light met- 
als, for which the approximation of free electrons is well 
justified, could be candidates for the very transparent 
interface between them. If aluminum is chosen as a su- 
perconductor (ei?(Al) = 11.7 eV [lb]) then the series of 
noble metals as a counter-electrode (from £i?(Ag) = 5.5 
eV to £i?(Cu) = 7 eV) corresponds to Rf — 0.4 — 0.6, 
which implies essentially reduced interface transparency. 
It seems that perfect interface transparency can not be 
realized in an experiment. In the view of the intrinsi- 
cally reduced transparency of the interface between any 
two metals we conclude that our quasiclassical boundary 
resistance can be used in the full range of experimentally 
attainable interface transparencies. On the other hand, 
all theories which essentially exploit perfect transparency 
of the interface have to be re-examined to check survival 
of predicted effects against the reduced transparency of 
the interface. 
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Figure Captions 

Fig. 1. Dependence of the normalized interface con- 
ductances on the Fermi energy ratio Rf- The scale of 
the ordinate is logarithmic. 
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